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We present high-resolution optical birefringence measurements upon sequential filling of an array 
of parallel-aligned nanochannels (14 nm mean diameter) with rod-like (acetonitrile) and disc-like 
(hexafiuorobenzene) molecules. We will demonstrate that such birefringence isotherms, when per- 
formed simultaneously with optically isotropic and index-matched counterparts (neopentane and 
hexafluoromethane) , allow one to characterize the orientational state of the confined liquids with 
a high accuracy as a function of pore filling. The pore condensates are almost bulk-like, optically 
isotropic liquids. For both anisotropic species we find, however, a weak orientational order (of a few 
percent at maximum) upon film-condensation in the monolithic mesoporous membrane. It occurs 
upon formation of the second and third adsorbed layer, only, and vanishes gradually upon onset of 
capillary condensation. Presumably, it originates in the breaking of the full rotational symmetry of 
the interaction potential at the cylindrical, free liquid-vapor interface in the film-condensed state 
rather than at the silica-liquid interface. This conclusion is corroborated by comparisons of our 
experimental results with molecular dynamics simulations reported in the literature. 

PACS numbers: 68.03.Fg, 68.08. Be, 05.70.Np 



I. INTRODUCTION 

Liquids confined in tiny pores a few nanometers across 
can have pr oper ties which differ markedly from their bulk 
counterpart^^!] Both the interaction with the confining 
walls and the bare geometrical confinement can induce al- 
tered molecular packings, layering, changes in the molec- 
ular dynamics as well as preferred collective orientational 
order in the case of liquids consisting of anisotropic build- 
ing blocks. 

Whereas both altered packing and layering- ^ as well 
as orientational ordei^^ has been established and sys- 
tematically probed for a sizeable number of molecular 
liquids confined in a semi-infinite manner at the free sur- 
face of liquids and solids, such order phenomena upon 
mesopore confinement are still elusive. In fact, except for 
the well documented cases of stroiigly anisotropic meso- 
genes, which for m n ematic phases^l, only a few experi- 
mental evidence^^ of such orientational order in meso- 
pores have been reported so far for anisotropic building 
blocks, that form solely isotropic bulk liquids. 

Due to the complex pore structure, tortuous pore ge- 
ometries and unknown pore wall/liquid interaction pa- 
rameters it is particularly challenging to examine orien- 
tational order in mesopores. Nowadays, however, the 
advent of tailorable mesoporous substrates with aligned, 
non-interconnected tubular pores allows one to study 
this phenomenology. Here we present a high-resolution 
optical birefringence study upon sequential filling of an 
array of parallel-aligned silica nanochannels via the va- 
por phase. We will demonstrate that such birefringence 
isotherms, when employed simultaneously for optically 



isotropic and anisotropic building blocks (with identical 
refractive index), allow one to characterize the orienta- 
tional state of the confined liquid with a high accuracy 
as a function of pore filling. 

We have chosen the disc-like fluorobenzene (CeHsF, 
referred hereafter as FB) combined with its optically 
isotropic counterpart tetrachloromethane (CCI4, TCM) 
and the rod-like acetonitrile (CH3CN, AN) in combi- 
nation with its optically isotropic partner neopentane 
(C5H12, NP), see Fig. [l] for illustrations of the molec- 
ular structures. The bulk liquids of FB and TCM have 
practically identical refractive indices, and so do AN and 
NP (see Table 1). 

The birefringence of these liquids have been measured 
as a function of fractional filling / of the nanochannel 
array over the full range from the empty matrix (/=0) 
to the completely filled case (/=!). For each pair the ex- 
perimental data on the anisotropic system are compared 
with the data of the isotropic system in order to extract 
a parameter typical of the collective orientational order 
of the molecules as a function of pore filling. 

In a former study we employed an identical technique 
to obtain information on confined liquid n-hexane. How- 
ever, we did not perform a comparison with an isotropic 
molecular counterpart, and hence had to use an effec- 
tive medium model in order to scrutinize for any collec- 
tive molecular orientation upon pore-condensation. It is 
understood that such a comparison can be affected by 
all kind of differences between the ideal model assump- 
tion and the real pore system, most prominently by a 
non-cylindrical pore shape and pore wall roughness. For 
n-hexane we did not find any clear cut evidence of a col- 
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FIG. 1: Volumetric sorption isotherms of the four molecular species invesitgated (as indicated in the figure), recorded at 
T =281 K in a monolithic mesoporous silica membrane. Open and solid symbols refer to ad- and desorption, respectively. The 
ray-tracing pictures illustrate the four molecular structures. The arrows represent the orientation of the principal axes of the 
polarizability tensors. Their lengths are proportional to the magnitude of the polarizabilities in those directions. 



lective orientational order. This will be contrasted by the 
systems studied here. 



II. EXPERIMENTAL 

Mesoporous Si has been prepared by electrochemical 
etching of a highly p-doped (100) wafer with a resistiv- 
ity of 0.01 il-cm by applying an electrical current of 12 
mA/cm^ for 8 hours. The resulting porous layer was re- 
leased from the bulk Si underneath, by increasing the 
etching current by a factor of ten. Subsequently the 
porous Si sheet was thermally oxidized under standard 
atmosphere. The oven temperature was ramped upwards 
from 100°C to 700°C over 36 hours, kept at 700 °C for 
48 hours and finally ramped downwards to room temper- 
ature over 12 hours. After this procedure the membrane 
was entirely oxidized. Standard analysis of a N2 sorp- 
tion isotherm, measured at a temperature of 77 K, gives 
a porosity of 41% and a pore diameter 2Rq of 14 nm. 
The sheet is 385 /im thick, hence the aspect ratio of the 
pores is ~ almost 2.5 -10^. Electron micrographs of the 
channels in the original porous silicon membranes^" in- 



dicate sizeable 1.0 ± 0.5 nm mean square deviations of 
their surfaces from an ideal cylindrical form. Presum- 
ably, the roughness in the oxidized samples are of similar 
magnitude. In order to minimize the influences of possi- 
ble pore geometry variations on our experimental results 
upon usage of different samples, the identical sample was 
employed throughout all experiments. 

Because of the perpendicular alignment of the chan- 
nels' long axes to the surface, the wafer is an effective 
medium that shows an uniaxial "shape" anisotropy, even 
in case the constituents have isotropic permittivities and 
refractive indices. Such examples of shape anisotropy 
have been observed in various systems, such as aligned 
carbon nanotube filmi^]^ porous Sp^'i^or semiconductor 
nanowireJi^. 

We have used a high-resolution optical polarization 
method for the accurate determination of the retarda- 
tion R at an incident angle of 33.7 deg. Details of the 
experimental setup can be found in Refs.^i^Hi^. The sam- 
ples have been kept in a closed optical chamber, held 
at a constant temperature (T = 281 K) and small vol- 
umetrically controlled portions of gas have been sub- 
sequently added (for adsorption) or removed (for des- 
orption) from the chamber by means of a gas handling 
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system equipped with a capacitive membrane gauge of 
100 mbar full scale. The equilibration time was one hour 
for each step. The pore condensate is specified by the re- 
duced pressure p = P/Pq and filling fraction / — N/Nq. 
Pq is the saturated vapor pressure of bulk hexane at given 
temperature T, N is the number of adsorbed molecules in 
the porous substrate and Nq is the number of molecules 
necessary for a complete filling of the pores. Optical bire- 
fringence measurements were only possible upon adsorp- 
tion in the hysteretic part of the isotherm. Upon des- 
orption a sizeable light scattering typical of a coarsening 
of the pore condensate upon desorption (see a discussion 
of this phenomenology in Ref.^i^ and references therein) 
hampered polarimetry measurements. The observed re- 
tardation R has been converted into the optical birefrin- 
gence An by means of the formalism described in detail 
in Refs.liSUei. 



III. RESULTS AND DISCUSSION 

A. Sorption Isotherms 

The volumetric sorption isotherms of all four molecu- 
lar species have a shape typical of condensation of van- 
der-Waals molecules in mesoporous substrates, see Fig. 
[T] The initial reversible part is due to adsorption of a 
film on the pore walls. The hysteretic part between ad- 
sorption and desorption is characteristic of capillary con- 
densation and evaporation. The subsequent plateau in- 
dicates the completion of pore filling, whereas the final 
increase is due to the condensation of vapor on the ma- 
chined rough metal walls of the sample cell and fina lly to 
bulk condensation. According to theoretical modelJ^^'^ 
for a cylindrical pore of uniform cross section the adsorp- 
tion branch should be vertical in the regime of capillary 
condensation. The finite slope of both branches is typi- 
cal of a broad distribution of pore diameters. Because of 
the rather gradual increase of the adsorption branch, the 
filling fraction fc at which capillary condensation starts 
cannot be unambigiuously determined. Judging from our 
experience with similar van-der-Waals systems condensed 
in more regular cylindrical pores, e.g. SBA-15, we will 
just refer to an approximate value of 0.5 for fc^ . 

B. Birefringence isotherms 

In Fig. [2] the birefringence An of all four molecular 
species is plotted as a function of filling fraction / as 
recorded upon liquid adsorption. The general shape of 
these curves (i.e. the finite starting value of An, the 
maximum at / about 0.3, and the subsequent decrease to 
very low values at / approaching 1) is similar to that ob- 
served previously for hexan^^^. It is consistent with what 
is expected from a model describing aligned pores filled 
with an isotropic medium in an optically isotropic ma- 
trix. Nevertheless there are slight differences 5 between 



the An(/) data of the optically anisotropic molecules and 
their isotropic counterparts, (5=Ananiso — Aniso. In Fig. 
[3] we show S as function of / for both the AN/NP and 
the FB/TCM pair. 

For the AN, the molecule with the prolate polarizibil- 
ity tensor, Atian > Atinp, that is S positive, whereas for 
FB with the oblate tensor, AnpB < AnxcM, <5 < 0. This 
is evidence of an intrinsic orientational order of the AN 
and FB pore liquids, beyond the optical "shape" bire- 
fringence brought forth by the parallel alignment of the 
pores. In the present context birefringence is sensitive 
to cos^ 6 where 6 is the angle between the pore axis and 
the direction of maximum (minimum) polarisibility a of 
the AN (FB) molecule (We ignore the small difference 
between the components and ai, of FB). Moreover, 
given the small and rigid character of the molecules em- 
ployed along with the weakness of any of the interfacial 
potentials when compared to the covalent intra-molecular 
bondings, we feel encouraged to believe that no addi- 
tional anisotropics within the single molecules are in- 
duced by the interfacial potentials, and thus hamper our 
analysis. 

An idea of the degree of the preferred orientation 
is obtained from the comparison with a calculation of 
the birefringence An for complete orientational order 
of these pore liquids, that is = (see Ref.^^ for de- 
tails of the model employed). In order to account for 
the observed signs of S, the AN molecules are all ori- 
ented with their axis parallel to the pore axis, whereas 
for FB it is the direction of smallest polarizibility that 
is parallel to the pore axis. For both cases the calcu- 
lated birefringence changes strongly and monotonously 
as function of /, starting at An =0.018 at / = and 
attaining a value of 0.17 at / = 1 for AN and -0.27 
for FB. Calling the difference between these calculated 
An- values and the experimental ones 5ma.x, the quantity 
q — S/Sniax — l/2(3cos^6'— 1) serves as orientational 
order parameter. The results on q{f) are included in 
Fig. [3} 

The values of q are small, the maximum value of q is 
of the order of one percent for both systems, and even 
these small values are only observed in a restricted /- 
range and by no means in the full range of /. Thus, 
ignoring a quite unlikely, peculiar angle of collective ori- 
entational order (discussed below) , we can conclude that 
the AN and FB pore liquids are almost bulk-like, opti- 
cally isotropic liquids. 



C. Discussion 

The residual preferred orientation can in principle re- 
sult from a breaking of the full rotational symmetry of 
the bulk liquid by the pore walls and/or the liquid- vapor 
interface. 

We concentrate on AN for which the structure and 
the dynamics of the liquid at the liquid-vapor interface 
and at a silica wall have been studied by molecular dy- 
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refractive indices 


molecule 


polarizability a 


vapor pressure po (8) 


71 


na,nb,nc 


^0 ,71e 


ecu - tetrachlormethane - TCM 




/7.474 0.0 0.0 \ 

0,0 7.474 0,0 
\ 0.0 0.0 7.474/ 




68.1 mbar 


1.459 


isotrope 




CeHsF - fluorobenzene - FB 




10.501 0.0 0.0 
0.0 10.744 0.0 
0.0 0.0 3.258 


) 


43.6 mbar 


1.463 


71^ = 1.584 
Tib = 1.601 
Tic = 1.161 


7io — 1.593 

Tie = 1.161 


C5H12 - neopentane - NP 




/8.013 0.0 0.0 \ 

0.0 8.013 0.0 
\ 0.0 0.0 8.013/ 




960.5 mbar 


1.342 


isotrope 




CH3CN - acetonitrile - AN 




/2.581 0.0 0.0 \ 

0.0 2.581 0.0 
\ 0.0 0.0 5.048/ 




51.8 mbar 


1.346 


Ua ^ nb = 1.248 
Uc = 1.524 


7io = 1.248 

Tie = 1.524 



TABLE I: Saturated vapor pressures (at T =281 K) and optical data sets of the molecules investigatecP^. 



namics simulation£p2H24|_ According to this work the AN 
molecules of the liquid next to the wall show a broad 
distribution of tilt angles 7 with respect to the normal 
but with almost vanishing probability for 7=0, thus the 
molecules are preferentially aligned perpendicular rather 
than parallel to the wall. In the present experiment such 
preferred orientations should lead to a negative contribu- 
tion to 5 that should develop already in the submonolayer 
region of the sorption isotherm and persist up to com- 
plete filling. By contrast, we find vanishing values of 5 in 
that regime. Hence, the experiment is compatible with 
an orientationally disordered state of the first molecules 
adsorbed. See the low starting values of q{f ), Fig. [sj that 
arise from the fact that the birefringence curves An{f ) of 
AN and its isotropic reference liquid start with the same 
slope (Fig. [2]). Of course pore walls of our samples are 
expected to be rough, both on a microscopic and on a 
mesoscopic scale. We anticipate that the resulting dis- 
tribution of the surface normal can easily wipe out any 
preferred orientation with respect to the macroscopically 
defined direction of the pore axis. Note that, however, 
also for e = 54.73 deg (Scos^ 6-1 = 0), the so-called 
magic angle, the resulting q and thus also S vanishes. 
Thus, in principle, we cannot exclude by our experiments 
values of 7 identical or close to the magic angle. Albeit, 
in the light of the aforementioned molecular dynamics 
studies this peculiar configuration is not expected--. 

In the following we argue that the finite, positive val- 
ues of the orientational order parameter observed in the 
experiment are induced by the free surface of the pore liq- 
uid. According to the MD simulations the AN molecules 
at the liquid/vapor- interface have a slight preference for 
parallel alignment (7=0). For cylindrical liquid/vapor- 
interfaces parallel to the pore axis this translates into 
positive values of 5, in agreement with the experiment. 
One also expects that such interfaces are much smoother 
than the pore walls, thereby introducing less uncertainty 



when changing over from ^ to 9. 

Finite values of q occur for intermediate values of / 
ranging from about / = 0.1 to / = 0.6. The delayed 
onset of q{f) at low / is related to the fact that the free 
surface of the liquid can only develop for liquid films on 
the walls that are sufficiently thick. According to Hu and 
Week^2llthe AN molecules within 2.5 nm from the silica 
wall are practically immobile. In the present samples 
with rough and highly curved pore walls the thickness, 
above which a free liquid/vapor surface can develop, is 
presumably even thicker. 

The vanishing values of g(/) beyond / = 0.7 is a more 
problematic point. In the regime of capillary condensa- 
tion (/ > /c = 0.5) the area of the cylindrical liquid- 
vapor interfaces is expected to approach zero at / = 1 
in a linear fashion. Thinking exclusively in such terms, 
q{f) should show the same /-dependence. It is, how- 
ever, well established that in the regime of capillary con- 
densation va por b ubbles coexist with completely filled 
pore sectiongiSSSl, Thus there are not only cylindrical 
but also meniscus-shaped liquid-vapor interfaces. With 
the molecules still preferentially parallel to the interface, 
the latter type interface supplies negative contributions 
to 6. The fact that the partially filled silica pores with 
/-values within the capillary condensation regime are op- 
tically tran sparen t suggests that the vapor bubbles are 
pretty smalpSESl^ but nevertheless it is surprising that 
the menisci are so abundant that their contribution to 
the preferred orientation of the molecules can compen- 
sate that of the cylindrical interfaces. Their great num- 
ber could originate in the sizeable pore size variation and 
the corresponding plethora of nucleation sites for liquid 
bridges within the channels. 

For FB finite values of q are observed in an /-range con- 
siderably narrower than for AN (Fig. [s]). The negative 
values of 6 are related to a preference of the normal of the 
aromatic ring pointing along the pore axis which trans- 
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FIG. 2: (Color online) The optical birefringence An vs / = N/Nq upon adsorption as determined from the measured phase 
retardation R for NP (solid circles), AN (solid squares) in the left panel and for FB (solid triangles) and TCM (solid circles) 
in the right panel. The insets show the low-/ data on an expanded An-scale. 



lates into a perpendicular orientation with respect to the 
pore walls and to the cylindrical liquid/vapor interfaces. 
As for AN the vanishing value of q{f) for f 's, / > 0.3, and 
for / approaching 0, suggests that the present experiment 
does not produce evidence for pore wall induced prefer- 
ential alignment, it rather suggests that any alignment 
observed is due to the liquid/vapor interface. Neverthe- 
less it is remarkable that the intrinsic birefringence is 
confined to such a narrow /-range which corresponds to 
roughly the second and third monolayer adsorbed on the 
pore walls. 

In liquid CgHg (and also in CeFg) the quadrupolar in- 
termolecular interactions favor configurations where the 
symmetry axes of neighbouring molecules arrange per- 
pendicular (f^Sl and references therein) . We assume that 
this interaction is also dominant in FB (in spite of its 
weak dipole moment) and that these configuration even 
prevail at the free liquid/vapor interface. In such a situ- 
ation the combined contribution of neighbors to q is just 
zero, thereby explaining why the FB molecules at the 
free surface do not contribute to the birefringence. We 
speculate that finite contributions can just occur in the 
special situation that such "T"-like pair configurations 
have not yet developed. This is likely to be the case at 
/ about 0.1, where a fractional monolayer exists on the 
immobile molecules next to the wall. 



Coasne et al. found in molecular dynamics simula- 
tions that the orientation of CeHg at the silica wall is 
very sensitive to its hydroxilation'^^. Whereas for highly 
hydroxylated pore walls the ring is arranged perpendic- 
ular to the surface, partially hydroxylation leads to a 
disordered state in the first monolayer. Infrared spec- 
troscopy measurements indicate that the pore walls are 
highly hyd roxyla ted after the preparation procedure em- 
ployed her^SlMl should expect wall-induced orienta- 
tional order, in contrast to our observation. Therefore, 
we have to resort to the assumption, similar to our con- 
siderations regarding the experiments on AN, that the 
orientational order in the wall proximity is washed out 
by roughness. 



IV. SUMMARY 

We presented a high-resolution optical birefringence 
study upon filling of a mesoporous monolithic silica mem- 
brane by four molecular species. The rod-like molecules 
acetonitrile and its optically isotropic and refrative- 
index-matched counterpartner neopentane as well as the 
disc-like fluorobenzene and its optically counterpart- 
ner tetrafluoromethane. A comparison of birefringence 
isotherms allowed us to gain detailed information on the 
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FIG. 3: 5 and q vs. / for the rod-like acetonitrile (squares) 
and the disc-hke fluorobenzene (triangles) condensed into a 
mesoporous silica membrane. 

collective order of the anisotropic building blocks, with- 
out having to resort to detailed model calculations on the 
changes of the optical birefringence upon capillary filling 
of the optically birefringent mesoporous matrix. 

We find a very weak orientational order. The pore- 
condensed liquid is bulk-like, optically isotropic. Pre- 



sumably, the weak orientational order is imposed by the 
breaking of the full rotational symmetry in the interac- 
tion potential at the cylindrical liquid/vapor interfaces 
upon film-condensation. Surpringly, upon onset of capil- 
lary condensation this preferred orientational order van- 
ishes quickly which can only partially be understood by 
the vanishing liquid/ vapor interfaces and the formation 
of concave menisci. 

In general, the interpretation of our results signifi- 
cantly profited from molecular dynamics simulations in 
planar and cylindrical pore geometry, available in the lit- 
erature. Thus our study also demonstrates that a thor- 
ough interpretation of birefringence isotherms crucially 
relies on comparisons with such calculations. Alterna- 
tively comparisons with experimental results concerning 
molecular orientational order at the free, planar liquid 
surface, gained for example by vibrational sum frequency 
spectroscopy^, will help to reduce the intrinsic ambigu- 
ities in the interpretation of the experimental method 
presented. 

For the future, we suggest to systematically investi- 
gate the influence of pore wall roughness and chemistry, 
e.g. by silanization of the silica walls, or the usage of al- 
ternative, optically transparent mesoporous substrates, 
for example monolithic alumina membranes. The result- 
ing variation of the pore wall/fluid interaction and re- 
duced geometric pore wall irregularities may allow one to 
find and explore also collective orientational order at the 
fiuid/solid interface, as expected by molecular dynamics 
studies. 
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